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ABSTRACT 

Time-dependent and direction-dependent neutrino and gravitational-wave (GW) signatures are presented for a set of three-dimensional 
(3D) hydrodynamic models of parametrized, neutrino-driven supernova explosions of non-rotating 15 and 20 Mq stars. We emploed 
an approximate treatment of neutrino transport based on a gray spectral description and a ray-by-ray treatment of multi-dimensional 
effects. Owing to the excision of the high-density core of the proto-neutron star (PNS) and the use of an axis-free (Yin- Yang) overset 
grid, the models can be followed from the post-bounce accretion phase through the onset of the explosion into more than one second 
of the early cooling evolution of the PNS without imposing any symmetry restrictions and covering a full sphere. Gravitational 
waves and neutrino emission exhibit the generic time-dependent features already known from 2D (axi-symmetric) models. Violent 
non-radial hydrodynamic mass motions in the accretion layer and their interaction with the outer layers of the proto-neutron star 
together with anisotropic neutrino emission give rise to a GW signal with an amplitude of ~ 5 - 20 cm in the frequency range of 
100-500 Hz. The GW emission from mass motions usually reaches a maximum before the explosion sets in. After the onset of the 
explosion the GW signal exhibits a low-frequency modulation, in some cases describing a quasi-monotonic growth, associated with 
the non-spherical expansion of the explosion shock wave and the large-scale anisotropy of the escaping neutrino flow. Variations 
of the mass-quadrupole moment caused by convective activity inside the nascent neutron star add a high-frequency component to 
the GW signal during the post-explosion phase. The GW signals exhibit strong variability between the two polarizations, different 
explosion simulations and different observer directions, and besides common basic features do not possess any template character 
The neutrino emission properties (fluxes and effective spectral temperatures) show fluctuations over the neutron star surface on spatial 
and temporal scales that reflect the different types of non-spherical mass motions in the supernova core, i.e., post-shock overturn flows 
and proto-neutron star convection. However, because very prominent, quasi-periodic sloshing motions of the shock caused by the 
standing accretion-shock instability are absent and the emission from different surface areas facing an observer adds up incoherently, 
the modulation amplitudes of the measurable neutrino luminosities and mean energies are significantly lower than predicted by 2D 
simulations. 

Key words. Supernovae: general — neutron stars: general 



1. Introduction 

The electromagnetic signature of core-collapse supernovae has 
been exploited comprehensively and thoroughly by countless 
observations during the past decades, providing only indirect in- 
formation about the explosion mechanism, however. The up to 
now only recorded neutrino signal of a core-collapse supernova 
(SN1987A) confirmed the idea that the collapse of the core of a 
massive star to neutron star densities provides the necessary en- 
ergy for the explosion (Baad e & Zwickvll934l) . Because gravita- 
' tional waves (GW), the only other means to probe the supernova 
engine besides neutrinos, are yet to be detected, supernova mod- 
elers are preparing for this prospective measurement by predict- 
ing the gravitational wave signature of core-collapse supernovae 
with ever increasing realism (for rev iews, see e.g., iKotake et aD 
120061: IOt3l20()9l:lFrver & Newll201 Ih . 

According the Einstein's theory of general relativity (GR), 
gravitational waves will be generated by any mass-energy distri- 
bution whose quadrupole (or higher) moment varies non-linearly 
with time. In core-collapse supernovae this criterion is satis- 
fied by time-dependent rotational flattening particularly during 
collapse and bounce, prompt post-shock convection, non-radial 
flow inside the proto-neutron star and in the neutrino-heated hot 
bubble, the activity of the standing accretion shock instability 



(SASI), asymmetric emission of neutrinos, and by asymmetries 
associate d with the effects of magnetic fields (for a recent review 
see, e.g., IOttjl2009l and references therein). While significant ro- 
tational deformation and dynamically relevant magnetic fields 
require particular progenitors that possess some (considerable) 
amount of rotational and magnetic energy or that must rotate 
fast and differentially (additional differential rotation develops 
during collapse) to amplify an initially weak magnetic field, all 
other processes are genuinely operative in any core-collapse su- 
pernova. 

Obviously, an adequate modeling of these effects and 
an accurate prediction of the GW signal ultimately requires 
three dimensional (3D) hydrodynamic simulations covering 
the collapse, bounce, and post-bounce evolution (at least) 
until a successful launch of the explosion and including a 
proper treatment of neutrino transport and the relevant micro- 
physics. However, most studies of the past thirty years were 
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iDimmelmeier et all 20081). or were restricted to axisymmetric 
(2D) models ('Muller et al."2004'; 'Ott et al."2006'; 'Kotake et al.' 
2007; Mareketal. 2009; Murphy et al. 2009; Yakuninetal. 
2010l) . Several authors also investigated the influence of mag- 



netic fields on the GW signal during the collapse and earl y post- 
bounce evolutiori assu r ning axisymmetry ( Kotake et al.l 120041; 
'Yamada & Sawai"2004 lKotake et al.ll2005 ; Oberaaulinaer et al.' 
2006 a. b) an d no symmetry restriction at all (Scheidegger et al. 



20081 I20l0h . The GW signal c aused by aspher ical neutrino 
emission was first studied by 'Epstein' ('l978h and subse- 



quer itly by [Bu rrows & Hayes (1996), Muller & Jaiikal (Il997l) . 
and iKotake eTaLl (2007, 2009a b, 2011), where the investiga- 
tions bylMMeE&Janka (.1997) and Kotake et al. (2009b. 2011.) 
considered also 3D, i.e., non-axisymmetric models. 

Concerning 3D post-bounce models, the topic of the study 
presented here, Muller & Janka ( 1997) were the first to analyze 
the GW signature of 3D non-radial flow and anisotropic neutrino 
emission from prompt post-bounce convection in the outer layers 
of a proto-neutron star during the first 30 msec after supernova- 
shock formation. Because of smaller convective structures and 
slower overturn velocities, the GW wave amplitudes of their 
3D models are more than a factor of 10 lower than those of 
the corresponding 2D ones, and the wave amplitudes associ- 
ated with anisotropic neutrino emission are a factor o f 10 higher 
than those cause d by non-radial gas flow. Fryer (2004) and 
iFrver et aP (l2004l) presented gravitational wave signals obtained 
from 3D core-collapse simulations that extend up to 150 msec 
past bounce and were performed with a Newtonian smoothed 
particle hydrodynamics code coupled to a three-flavor gray flux- 
limited diffusion neutrino transport scheme. Gravitational wave 
emission occurs owing to matter asymmetries that arise from 
perturbations caused by precollapse convection, core rotation, 
and low-mode convection in the explosion engine itself, and ow- 
ing to anisotropic neutrino emission. Kotake et al. (2009b) sim- 
ulated 3D mock-up models that mimic neutrino-driven explo- 
sions aided by the SASI, and computed the GW signal resulting 
from anisotropic neutrino emission by means of a ray-tracing 
method in a post-processing step. They pointed out that the grav- 
itational waveforms of 3D models vary much more stochas- 
tically than those of axisymmetric ones, i.e., in 3D the GW 
signals do not possess any t emplate character. H owever, when 
considering rotating models, iKotake et al.l (1201 ll) argue that the 
GW waveforms resulting from neutrino emission exhibit a com- 
mon feature, which results from an excess of neutrino emission 
along the spin ax i s due to the growth of spiral SASI modes. 
IScheidegger et al.l (l2008h simulated the collapse of two rotating 
and magnetized cores in 3D until several 10 msec past bounce, 
applying a parametrized deleptonization scheme ('Liebe ndorfeii 
[2005) that is a goo d approximatio n until a few milliseconds past 
bounce. [Scheideg ger et al.l (1201 Ol) extended their study by sys- 
tematically investigating the effects of the equation of state, the 
initial rotation rate, and both the initial magnetic field strength 
and configuration on the GW signature. They also simulated 
three representative models until ~ 200 msec of post-bounce ac- 
cretion (no development of explosions) incorporating a treat- 
ment for neutrino transport based on a partial implement ation of 
the is otropic diffusion source approximation (Liebendorf er et af] 
120091) . 

In the following we present the GW analysis of a set of 
parametrized 3D models of neutrino-powered supernova ex- 
plosions covering the evolution from core bounce until ~ 1 .4 s 
later, where the high-density inner core of the proto-neutron star 
(PNS) is excised and replaced by a time-dependent boundary 
condition and a central point mass. The neutrino transport is 



treated by an approximate solver b ased on a ray-by-ray treatment 
of the multi-dimensional effects (IScheck et al.l 12006 '). Because 
we analyze the GW signal arising from both non-radial mass 
flow and anisotropic neutrino emission, we discuss the neutrino 
emission of these 3D models as well, and particularly address 
its multidimensional properti es, some of which have previou sly 
been investigated in 2D by Janka & Monchmeyer (1989a oh, 
Ott et al. (20 08), Kotak e et al. (200_9a ), Marek & Janka (2009^, 
iMarek et all ( l2009h . and lBrandt et alj JSoTTh . 

Based o n 2D si mulations of rotational core-collapse, 

iJank a & Monc hme ved fl 9 8 9 a b ) pointed out that the energy out- 
put in neutrinos seen by an observer may vary as much as a 
factor of 3 with his inclination angle relative to the rotation 
axis, while for the neut rino energy much sm all er angular varia- 
tions are to be expected. iMarek et al.l ( l2009l) and lMarek & Janka[ 
(2009) found that neutrino luminosities and mean energies show 
quasi-periodic time variability in their 2D simulations of non- 
rotating and slowly rotating 15 Mq progenitors covering up to 
~700ms past bounce. Caused by the expansion and contraction 
of the shock in the course of SASI oscillations, the level of the 
fluctuations is $50% for the luminosities and roughly 1 MeV 
for the mean neutrino energies. The luminosity fluctuations are 
somewhat bigger for Ve and Ve than for heavy-lepton neutrinos. 
The neutrino quantities also depend on polar angle as a conse- 
quence of the preference of the SASI motions along the symme- 
try axis of the 2D models. Additional short-wavelength spatial 
variations of the average radiated energies and of the neutrino 
fluxes are caused by local downdr afts of ac creted matter. This is 
in accordance with the findings o f iMuller & Janka (1997), who 
inferred from a post-processing analysis of the neutrino emis- 
sion anisotropy that features in the GW signal of their 2D mod- 
els of convection in the hot-bubble region are well-correlated 
with structures in the neutrino signal. The features are associ- 
ated with sinking and rising lumps of matter and with temporal 
variations of aspherical accretion flows toward the proto-neutron 
star. Kotake et al.. (,2009a) computed neutrino anisotropics with 
a ray-tracing scheme by post-processing their 2D SASI models 
and derived analytical expressions for evaluating GW signals for 
neutrinos in 3D models, too. A generalization of these expres- 
sions will be presented below. iBrandt et al.l (1201 II) performed 
2D multi-group, multi-angle neutrino transport simulations for 
both a non-rotating and rapidly rotating 20 Mq model extending 
~400ms beyond bounce. Their simulations predict that the neu- 
trino radiation fields vary much less with angle than the matter 
quantities in the region of net neutrino heating because most neu- 
trinos are emitted from deeper radiative regions and because the 
neutrino energy density combines the specific intensities as inte- 
grals over sources at many angles and depths. The rapidly rotat- 
ing model exhibits strong, flavor-dependent asymmetries in both 
peak neutrino flux and light curves, the peak flux an d decline rate 
having pole-equator ratios :^3 and ^2, respectively. ^ Brandt et al.l 
(1201 ih also provide estimates of the detectability of neutrino 
fluctuations in IceCube and Super-Kamiokande as previously 
done by Lund et al. (2010,.) on the basis of the .Marek et al., (.20091) 
non-rotating models. 

The paper is organized as follows: in Section 2 we discuss 
the numerical methods, the input physics, and the properties of 
the progenitor models and the set of 3D simulations that we 
analyzed. Section 3 contains a description of the formalism we 
used to extract the observable neutrino signal of our 3D models, 
and a discussion of some of its properties relevant for the cor- 
responding GW signal. In Section 4 we give the formalism nec- 
essary to calculate the GW signature of 3D non-radial flow and 
anisotropic neutrino emission, and discuss the GW signature of 
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the investigated 3D models. Finally, in Section 5 we summarize 
our results and discuss shortcomings and possible implications 
of our study. 



Table 1. Some properties of the analyzed 3D models including 
the mass of the progenitor star Mps, the mass of the neutron star 
Mns, the time of explosion fexp, and the explosion energy fiexp at 
the time ff after bounce when we stopped the simulation. Note 
that 1 B = 1 bethe = 10^' erg. 



2. Model setup 

2.1. Code and computational grid 

The 3D supernova models we analyzed for their neutrino and 
GW signature have been simulated with the explicit finite- 
volume, Eulerian, multi-fluid hy drodynamics code Prometheus 
dFrvxeU et al]|1991l: iMuUer et al.Ll991a.b,) . This code integrates 
the multidimensional hyd rodyna r nic eq uations using the dimen- 
sional splitting method of Strang ('1968"), the piecewise parabolic 
method of Collela & Woodward ( 1984), and a Riemann solver 
for real gases proposed by iColella & Glad (Il985h . Inside grid 
cells with strong grid-aligned shocks fluxes computed fro m the 
Riema nn solver are replaced by the AUSM-H fl uxes o f Lioul 
(11996 *) in order to prevent odd-even decoupling (Ouirld fT994 ). 
Nuclear species are advect ed using the consistent multi-fluid ad- 
vection (CM A) scheme of lPlewa & Mulled d 19991) . 

The simulation code employs an axis-free overlapping "Yin- 
Yang" grid ( Kagevama & Sato. .2004.) in spherical polar coor- 
dinates, which was re cently implemented into Prom etheus, for 
spatial discretization ("Wongwat hanarat et al.ll2010al) . The Yin- 
Yang grid relaxes the CFL-timestep condition and avoids nu- 
merical artifacts near the polar axis. Concerning the ray-by-ray 
neutrino transport no special procedure needs to be applied for 
the Yin- Yang grid. The (scalar) quantities involved in the trans- 
port algorithm are computed on both the Yin and Yang grid and 
are then linearly interpolated as any other scalar quantity. 

The grid consists of 400(r) x 47(0) x 137(0) x 2 cefls cor- 
responding to an angular resolution of 2° and covers the full An 
solid angle. The radial grid has an equidistant spacing of 0.3 km 
from the inner grid boundary out to r = 80 km (models W15 
and N20; see Table 1) or 115 km (model L15; see Table 1), re- 
spectively. Beyond this radius the radial grid is logarithmically 
spaced. The outer radial grid boundary 7?ob is at 18000 km, which 
is sufficiently far out to prevent the supernova shock from leav- 
ing the computational domain during the simulated epoch. This 
radial resolution suffices the requirement that there are always 
more than 15 radial zones per decade in density. 

A central region, the dense inner core of the proto-neutron 
star (PNS) at p ^ 2q12...i3 gem"-', is excised from the compu- 
tational domain and replaced by an inner time-dependent radial 
boundary condition and a point mass at the coordinate origin. 
The radius of the inner radial boundary shrinks according to 
Eq.(13) of lScheck etal] (|2008|) . For the W15 and N20 models 
the initial and final (asymptotic) boundary radii are - 65 km 
and R\t^ - 15 km, respectively. For the L15 models the corre- 
sponding radii are 82 km and 25 km, accounting for a less ex- 
treme contraction of the neutron star within the simulation time. 
The timescale for the contraction is fib = 1 sec for all models. 
This choice of parameters implies Ri\, a: 19 km at 1.3 s for mod- 
els W15 and N20, and R^b ^ 30 km at 1.4 s for the L15 mod- 
els. Hydrostatic equilibrium is assumed at the inner radial grid 
boundary R\\„ which is thus a Lagrangian (co-moving) position, 
while a free outflow boundary condition is em ployed at the outer 
radial grid boundary (for rnore d etails, see I Won g wathanarati 
1201 lHWongwaflianarat et alJl2010ah . 



Model 


Mps 


Mns 


^exp 


F 

^exp 


U 


[Mo] 


[Mo] 


[ms] 


[B] 


[s] 


W15-2 


15 


1.37 


248 


1.13 


1.3 


W15-4 


15 


1.38 


272 


0.94 


1.3 


L15-2 


15 


1.51 


381 


1.74 


1.4 


L15-3 


15 


1.62 


477 


0.84 


1.4 


N20-2 


20 


1.28 


265 


3.12 


1.3 



2.2. Input pliysics 

Self-gravity is fully taken into account by solving Poisson's 
equation in i ntegral form using an expa nsion into spherical har- 
monics as in lMiiUer & Steinmetzl(ll995h . The monopole term of 
the potential is corrected for general relativistic effects as de- 
scribed in Schecketal. (2006) and Arcones et al. (2007). The 
cooling of the PNS is prescribed by neutrino emission properties 
(luminosities and mean neutrino energies) as bound ary condition 
at the inner radial grid boundary (for details, see IScheck et al.l 
'2006'). The contraction of the PNS is mimicked by the move- 
ment of the inner radial grid boundary (see Sect. l2.1b . "Ray-by- 
ray" neutrino tra nsport and neutrino -matter interactions are ap- 
proximated as in lScheck et al.l (l2006h by radial integration of the 
one-dimensional (spherical), gray transport equation for all an- 
gular grid directions (6, (p) independently. This approach allows 
for angular variations of the neutrino fluxes. The tabulated equa- 
tion of state (EoS) of iJanka & Miiller ( 1996) is used to describe 
the stellar fluid. It includes arbitrarily degenerate and arbitrarily 
relativistic electrons and positrons, photons, and four predefined 
nuclear species (n, p, a, and a representative Fe-group nucleus) 
in nuclear statistical equilibrium. 

2.3. Models 



We have analyzed a set of 3D simulatio ns (I Won g wathanarati 
1201 1[ IWongwaflianarat et ani2010bi l201 ih based on two 15Mo 
progenitor models (W15 and LI 5), and a 20 progen- 
itor model (N20). The W15 model is obtained from the 



non-rotating I5M0 progenitor sl5s7b2 of Wooslev & Weaver 



(1995), the L15 model from a star evolved by Limon gi et al 



(2000), and flie N20 model from a SN 1987A progenitor of 
Shigeva ma & Nomotol (Il990l) . The progenitor models were 
evolved through collapse to 15 ms after bounce with the 
Prometheus- Vertex code in one dimension (A. Marek and 
R. Buras, private communication) providing the initial models 
for the 3D simulations. To break the spherical symmetry of the 
initial models, random seed perturbations of 0. 1 % amplitude are 
imposed on the radial velocity (v,) field at the beginning of the 
3D simulations. Explosions are initiated by neutrino heating at 
a rate that depends on suitable values of the neutrino luminosi- 
ties imposed at the lower boundary chosen such that the desired 
value of the explosion energy is obtained. The evolution is fol- 
lowed until 1.3 s after bounce for the W15 and N20 progenitor 
models, while the L15 models are simulated until 1.4 s post- 
bounce. The GW analysis presented below comprises five mod- 
els (see Table 1), where models W15-2 and W15-4 differ only by 
the initial seed perturbations. The explosion energies, fiexp, given 
in Table 1 are instantaneous values at the end of the simulations 
(f = ff), adding up the total energies (kinetic + internal + gravi- 
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tational) in all zones where the sum of these energies is positive. 
The explosion time, fgxp, is defined as the time when this sum 
exceeds a value of 10'*'^ erg, roughly corresponding to the time 
when the average shock rad ius is 400 to 500 km (see, however, 
iPeicha & Thompsoiil (1201 lb for an alternative definition of the 
time of the onset of the explosion). 

3. Neutrino signal 

The non-radial motions caused by the SASI and convection in 
the neutrino-heated hot-bubble as well as by convection inside 
the proto-neutron star (driven by Ledoux unstable lepton gradi- 
ents) give rise to a time-dependent, anisotropic emission of neu- 
trinos of all fla vors, and thus to the emission of gravitational 
wav es (Epstein 1978; Burrows & Haves 1996; Mii ller & Jankal 
[1997: Kotake et al. 2007. 2009a.b) ." as discussed in Sect.g] 
We have analyzed this emission for the 3D models discussed 
above (see Sect. l2.3l l. particularly addressing its multidimen- 
sional properties. 




Fig. 1. Sketch illustrating the various quantities involved when 
deriving the observable luminosity of a radiating source, whose 
visible surface is shaded in blue. When the observer is at infinity, 
the emitting surface is a sphere, and is measured from the 
origin of this sphere, one obtains a = y. 



3.1. Formalism 

To derive observable luminosities of an emitting source we con- 
sider an observer located at far distance D from that source (see 
Fig lD- According to definition the flux measured by the ob- 
server is given by the following integral at the position of the 
observer: 

FiDj) ^ i^dcjiiIiD,co,t), (1) 

where fi is the cosine of the angle between the direction of the ra- 
diation and the line of sight (between the observer and the center 
of the source), a> denotes the radiation direction at the observer's 
location (defined by a pair of angles), and doj is the solid angle 
element around the radiation direction ai. The intensity / adopts 
non-zero values within the opening angle subtended by the emit- 
ting surface (not necessarily a sphere). We note that here and in 
the following we suppress the dependence of the intensity on the 
neutrino energy and assume that energy-integrated quantities are 
considered (the outlined formalism, however, is valid also for an 
energy-dependent treatment). The integration over ai at the ob- 
server's location can be substituted by an integration over the 
emitting surface of the source, because the radiation intensity is 
constant along rays, i.e., 

I{D,oj,t)^Io(Ro,oJo,t) (2) 

for any ray arriving at the observer from the source (and zero 
otherwise), where /?„ denotes the position of a surface element 
of the emitting surface in the coordinate frame of the source and 
lOo the direction of the radiation field at that position toward the 
observer Note that we ignore in Eq. (|2]) the trivial effect that the 
time t for /„ relative to the time for / is subject to a retardation. 

Moreover, in the following we disregard spectral and angular 
corrections that may be relevant when the emitting surface is in 
relative motion to the observer or sitting deep in the gravitational 
potential of a compact star (in which case general relativistic 
(GR) effects like redshifting and ray bending would be impor- 
tant). Considering the source to be at rest is a good assumption 
for the neutrinospheric region in the supernova core after bounce 
(the velocities of mass motions in this layer are unlikely to be 
higher than some 1000 km/s, i.e., at most around one percent of 
the speed of light), while GR energy redshift is certainly of rele- 
vance on the ~ 10-20% level during the proto-neutron star cool- 
ing phase (f > 1 s after bounce), but much lower in the accretion 



and shock revival phases, when the forming neutron star is still 
considerably less compact than the final remnant. 

For a distant observer D » max{|/?„|) holds, i.e., the value of 
jj. is very close to one for the whole emitting surface. Denoting 
the solid angle subtended by a surface element of the emitting 
surface by dw, we have dw = dA^/D^, where dA^ = cosydA 
is the projected area of this surface element perpendicular to the 
line of sight, when y is the angle between the normal unit vector 
Ha of the emitting surface element dA and the line of sight (see 
Fig.[TJ, but specified to the case D » \Ro\. Hence, we obtain for 
the observable luminosity the expression 

Lo(f) = 47tD^F{D, t) = 4n I dA cos jIo{Ro, co^, t) . (3) 

«7vis.suif . 

In order to evaluate the integral in Eq. one needs to 
know the intensity !„ as a function of energy, emission direction, 
and time at every point of the radiating surface of the source. 
Determining I„{Ro, oi,,, in general requires calculating full- 
scale neutrino transport. With this quantity as the solution of the 
transport problem at hand, Eq. (|3]l can be evaluated directly (in 
general with energy dependence) by performing the integration 
over any surface that encloses the radiation (neutrino) emitting 
source and that lies outside the volume where radiation interacts 
with matter {i.e., the intensity /„ in all points on the chosen sur- 
face must be given in the reference frame of the observer and 
must be preserved on the way from the emission point to the 
observer). 

Our ray-by-ray transport approximation, however, yields 
only the local neutrino energy density EiR,,, t) and the neutrino 
flux density FiR,,, t). To estimate the neutrinos radiated from ev- 
ery point of the neutrino sphere to the observer, we therefore have 
to develop an approach that yields a reasonable representation of 
the direction-dependent intensity as function of the quantities de- 
livered by our transport approximation^. To this end, we assume 
that the neutrino distribution is axisymmetric around the normal 

' Note that the procedure described in the following does not depend 
on whether the transport is performed in the gray approximation or 
is energy dependent. We therefore suppress the energy variable in all 
transport quantities and introduce mean energies in our gray treatment 
instead of considering neutrino energies as a function of spectral fre- 
quencies. 
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vector Ha at all points Ro- This implies that the direction depen- 
dence of the intensity lo is described by the direction angle y 
only (see Fig.[Tli, i.e., !„ - hiRo, 7, t), and that the flux direction 
is given by «a- Assuming further that loij) can be approximated 
by the lowest two terms of an expansion in spherical harmonics 
(as in the diffusion approximation), one can write 



/o(/?o,r,f) = 



Fo(Ro,t) 
2n 



1 + - cosy 



(4) 



Because the radiation flux density is defined as the first an- 
gular moment of the intensity, one can easily verify that the nu- 
merical coefficient 3/2 of the dipole term ensures that the flux 
density Fo{Ro, t) (normal to the emitting surface element dA) is 
given by ^~ dtp d cos y cos yIo(Ro, J, t), if the radiating sur- 
face does not receive any incident neutrinos from outside (i.e., 
I„{R„, y, t) = for cos y < 0)3 

Inserting Eq. (|4| into Eq. (O, we find for the observable neu- 
trino luminosity the expression 



Lo(f) = 2 ( dAcosrF„(/?„,r) l + -cosr . (5) 

vis. surf. \ / 

We further define an observable mean neutrino energy according 
to 

<£)„(0 = ^, (6) 



where 



«7 VI 



Ln,o{t) 



Fo(R,„t) / 3 
dA cosy 1h — cosy) (7) 

VIS. surf . 



is the observable neutrino number flux with e being the mean 
energy of the neutrino energy spectrum radiated from point R,,. 

Our 3D radiation hydrodynamics code computes the time- 
dependent neutrino energy flux density, Fo{_Ro,t), and neutrino 
number flux density, F„^o(Ro, t), through a sphere of radius Ro = 
\Ro\ in dependence of the angular position £1 = (9, (p), but actu- 
ally stores the related quantities 



A(n,t) = 4jiRlF„(Ro,n,t), 



and 



A„(n,t) = 



A(£l, t) 
e(£l, f) 



AnRlF„JRo,n,t), 



(8) 



(9) 



because these quantities are constant in the ray-by-ray approxi- 
mation of the free streaming region. In this approximation, both 
the neutrino energy flux and the neutrino number flux are purely 
radial. 



^ Note that the requirement /„ > implies that Eq. l|4) is valid in 
the whole range of cosy > -|, which includes inward going radi- 
ation for cosy < 0. Extending the integration over all directions of 

validity, one obtams |^Fo = d^ J d cos y cos y/^ x Fg and 
Eo = d^ J d cos y /„ = ^Fo/c a; 3Fa/c. This means that we 

have « ^cEg, which is a reasonably good approximation of the re- 
lation between flux density and energy density at the neutrinosphere, 
where one typically obtains /„ = F„/(E„c) a | . . . 5 for the flux fac- 
tor at an optical depth between unity and about 2/3 in detailed neutrino 
transport calculations in spherical symmetry (Janka 1991). The ansatz 
of Eq. is therefore consistent with basic properties of the neutri- 
nospheric emission. Moreover, we note that the expression corresponds 
to the limb-darkening law /£(cosy)//£(l) = (2/5)(l -1- 3/2 x cosy) that 
can be derived on grounds of the Eddington approxim ation (see, e.g., , 
lMihalasl ( fT973) . page 61, or lMorse & FeshbacH d 19531) . page 187). 



Using Eqs. © and (|9]l, and the fact that dA - RldQ. with 
dQ = sin 9A6d(p for the special case of an emitting sphere of ra- 
dius R„, we can rewrite the general expression for the observable 
neutrino luminosity given in Eq. (|5]l in the form 



Loit) 



- r 

27r Jvis 



dQ cosyA(n,f) 1 1 + -cosyl (10) 



and that of the observable neutrino number flux given in Eq. (|7]) 
in the form 



L„,„{t) 



- r 

27r Jvi 



dO cosy A„(n,f) 1 + -cosy) , (11) 

vis. hem. \ 2 



where in both cases the integration is performed over the visible 
hemisphere. 

For the evaluation of the gravitational wave amplitude in 
Sect. l4.1.2l we will also need the quantity 



dA 
dQ 



(n,t) = Fo{RoM,t)Ri 



(12) 



and the corresponding angle-integrated quantity 
A(f) 



= (£ da^{a,t)^^S) dQA(n,f). (13) 

Jsurf . dQ 47r Xrf. 

For the later discussion of the results we finally define the 

surface-averaged neutrino flux density 

1 r 1 d£(r) 

{Fo){t)= — (h dnFo(Ro,nj)=-r-^-^, (14) 
4;r Jsurf. 47r7?2 df 



where 



d£(0 
df 



surf. 



dQ^(Q,f) = A(f) 
dQ 



(15) 



is the total energy loss rate at time f from the supernova core 
to all directions, which (because of the flux variations over the 
sphere) is no directly observable quantity. 

We have also analyzed the evolution of the neutrino flux 
asymmetry by calculating the angular pseudo-power spectrum 
of the neutrino energy flux variation 



A(£l, r) - A(f) 

m 



(16) 



where A{£1, t) and A(f) are defined in Eqs. ([8]l and ( fT3] ). respec- 
tively. The pseudo-power spectrum is given by the decomposi- 
tion of 5a(Q, f) in spherical harmonic coefficients 



Almif) 



dQ5A(Q,r)F,;(Q) 



(17) 



where J',*„(Q) is the respective (complex conjugate) spherical 
harmonics. For our mode analysis we actually used the pseudo- 
power coefficients Co = | AqoP and 



Ci = 



1 



21+ 1 



|A,( 



2£|A„„|- 



(18) 



for / > 0, respectively. 



^ Using of the ray-by-ray approximation has the advantage that the 
evaluation of the integrals on the rhs of Eqs.[5]and|7]does not require the 
specification of a suitable surface, but can be done on any sphere outside 
the neutrino-decoupling region, as Eqs.llOlandllllare independent of 
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Fig. 2. Snapshots of models W15-4 (left) and L15-3 (right) illustrating the four phases characterizing the evolution of our 3D 
models (see text for details). Each snapshot shows two surfaces of constant entropy marking the position of the shock wave (gray) 
and depicting the growth of non-radial structures (greenish). The time and linear scale are indicated for each snapshot. 
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quasi-spherical shock expansion phase 
shocl< formation and shocl< expansion 

pre-explosion phase 

growth of post-shocl< convection 
and of the SASI 

post-explosion accretion phase 
shock revival, shock acceleration, 
ongoing mass accretion 



post-accretion phase 
end of mass accretion, 
onset of a nearly spherical 
neutrino-driven wind 
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Fig. 3. Shock radius (top) and total {i.e., summed over all flavors) 
energy loss rate due to neutrinos (bottom) as functions of time 
for model W15-4. In the upper panel, the black curve shows the 
angle-averaged mean shock radius, the blue (red) curve gives 
the maximum (minimum) shock radius, and the vertical dashed 
line marks the time of the onset of the explosion as defined in 
Sect. l2.3l In the lower panel, the blue and red curves show the 
time evolution of An,ax(^^, and Amin(Q, t), the maximum and 
minimum value of A(Q, f) (Eq.[8]l on a sphere of 500 km radius, 
respectively. The black line gives the corresponding surface- 
averaged value A(f) (Eq.[T3]l. Note that the luminosities imposed 
at the inner radial grid boundary are kept constant during the first 
second and later are assumed to decay like r^^^. 



3.2. Resu/fs 
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The evolution of our models can be divided into four distinct Fig- 4. Same as Fig.[3]but for models L15-3 (uppermost two pan- 
phases (Figs.|2l[3]i. els) and N20-2 (lowermost two panels), respectively. 



7 



E. Miiller et al.: Parametrized 3D models 

(1) The first phase, the quasi-spherical shock-expansion phase 
(Fig. 121 top row), lasts from shock formation shortly after 
core bounce to 80 - 150 msec, when convection sets in. 
During this phase the shock rapidly propagates out to a ra- 
dius of ~ 200 km, where its expansion comes to a halt. 

(2) The second phase, the hydrodynamically vigorous pre- 
explosion phase, comprises the growth of post-shock con- 
vection and of the standing accretion shock instability, SASI 
(Fig.m second row from top). 

(3) The post-explosion accretion phase begins when energy de- 
position by v-heating in the post-shock layers becomes suf- 
ficiently strong to launch the explosion, and the total energy 
in the post-shock region ultimately becomes positive (see 
Sect. l2.3l for a definition). During this phase the shock accel- 
erates outward while gas is still accreted onto the PNS. This 
process is commonly called "shock revival" (Fig.|2] third row 
from top). 

Non-radial instabilities during the latter two stages cause 
considerable temporal and angular fluctuations of the neu- 
trino energy flux density as illustrated in Figs.|3]-|5] Besides 
the evolution of the shock radius, the figures show the 
surface-averaged neutrino light curve A(f), i.e., the total en- 
ergy loss due to neutrinos versus time (Eqs. [T3lfT5] ), together 
with the time evolution of the maximum and minimum val- 
ues of A(n, t) (Eq. ([8]i; the numerical evaluation is performed 
on an arbitrarily chosen sphere of 500 km radius). Distinct 
and high-amplitude spikes in A^n^iQ., t) are visible for sev- 
eral 100 msec and reflect violent post-shock convection, pos- 
sible SASI activity, and anisotropic accretion fluctuations af- 
ter the onset of the explosion. We have marked the explosion 
time fexp (see Section l231 and Table 1) by a vertical dashed 
line in Figs.|3]and|4] The post-explosion accretion phase lasts 
until -500 msec (models W15-4 and N20-2) or -700 msec 
(model L15-3) depending on the progenitor. 

(4) During the post-accretion phase, the fourth and final phase 
characterizing the evolution of our models (Fig.|2] bottom 
row), gas infall to the proto-neutron star has come to an end 
and the newly formed neutron star looses mass at a low rate 
in a nearly spherical neutrino-driven wind. We find consid- 
erably less temporal variability and a lower level of angular 
variation 10%) of the neutrino emission during this fourth 
phase (Figs.[3]-|5]l. 

While in model L15-3 the amplitudes of the neutrino emis- 
sion fluctuations decrease continuously, the other two mod- 
els exhibit growing temporal emission variations (though at 
a lower level than the earlier variability) during a later stage 
(notice the decrease/increase in A^ax- Ami,, in Figs.|3]and|4]l, 
which might be considered as a fifth evolutionary phase. This 
phase is associated with growing convective activity below 
the neutrinosphere. This PNS convection develops more or 
less strongly in the different models depending on the loca- 
tion of the convectively unstable region relative to the inner 
radial boundary of our computational domain. 

We have evaluated the time evolution of the neutrino energy 
flux asymmetry by producing 47r-maps that show the relative an- 
gular variation AFo/ {Fg) of the total (i.e., sum of all neutrino 
flavors) neutrino energy flux density across a sphere (normalized 
to the surface-averaged flux density; Eq.fT4li. Several snapshots 
of this evolution are shown for model W15-4 in Fig.|5] The evo- 
lution of the typical angular scales of the fluctuations is reflected 
by the pseudo-power spectrogram of the electron neutrino en- 
ergy flux variation (Eq.fTSTi in Fig.|6] top panels, which give the 
color-coded pseudo-power coefficient distribution normalized to 



neutrino-driven supernova explosions 




Fig. 5. Neutrino flux asymmetry at 170 msec, 200 msec, 
342 msec, 600 msec, and 1.3 sec (from top to bottom), re- 
spectively. The 47r-maps show the relative angular variation 
AFo/ {Fo) of the total {i.e., sum of all neutrino flavors) neutrino 
energy flux density over a sphere (normalized to its angular aver- 
age) for model W15-4. The maximum value is given in the lower 
right corner of each panel. Regions of higher emission are shown 
in bright yellow, while orange, red, green, and blue colors indi- 
cate successively less emission. Note that the color scale of each 
panel is adjusted to the maximum and minimum values at the 
coiTesponding time. The total energy loss rate due to neutrinos 
is given in the lower left corner. 
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Fig. 6. Pseudo-power spectrogram of the electron-neutrino energy flux density (top row) for models W15-4 (left), L15-3 (middle), 
and N20-2 (right), respectively. The panels in the middle row show the corresponding maximum pseudo-power coefficient C™" as a 
function of time, and the panels in the lower row give the relative angular variation of the electron neutrino flux density (maximum 
minus minimum flux density on the sphere divided by the angle-averaged flux density in percent) with time. 



the maximum value versus time. The variation of the pseudo- 
power coefficients with angular mode number is shown in Fig.|7] 
at selected times of 200 ms (blue), 400 ms (red), and 1000 ms 
(black). 

During the quasi-spherical shock expansion phase the level 
of angular fluctuations of F is low 10"^), while the fluctua- 
tion amplitudes of the total neutrino energy flux density reach a 
level of several 10% during the hydrodynamically vigorous sec- 
ond phase and the post-explosion accretion phase, where a few 
distinct regions or even single spots with an angular size of 10° 
to 20° dominate the emission (Fig.|5] panels 2 and 3). The mode 
number I of the dominant angular perturbation scale is of no rel- 
evance during the first phase, as the maximum pseudo-power co- 
efficient C™" (see Eq.[T8]l is tiny ^ 10"^ (Fig.|5] middle panels), 
i.e., the dominating I - 2 and / - 4 modes visible in the upper 
panels of Fig.|6]only reflect tiny angular perturbations imprinted 
presumably by the computational grid. When neutrino heating 
eventually causes significant non-radial flow during the second 
and third phases, C™" rises sharply to a level of ~ 10"^ (Fig.|6] 
middle panels), and the relative angular variations of the elec- 
tron neutrino flux density grow to the several ten percent level 
(Fig. 13 bottom panels). The latter quantity gives the maximum 
minus the minimum flux density on the sphere divided by the 
angle-averaged flux density in percent. Compared to the total 
neutrino emission in Figs.|3]-|5] the temporal and angular vari- 
ations in different directions are even more pronounced when 
considering the energy flux of the electron neutrinos or electron 
anti-neutrinos alone (Fig.|6]l, where angular variations can ex- 
ceed 100% in all models during the pre-explosion and accretion 
phases, and peak values are close to 200% during short episodes 
(Fig.|6] lower panels). 



During the vigorous pre-explosion phase including the post- 
explosion accretion stage, electron neutrinos and antineutrinos 
dominate the angular flux variations, while muon and tau neutri- 
nos (accounting for roughly 50% of the total luminosity) exhibit 
essentially isotropic emission in all directions. The reason of this 
finding is that and Vg are produced almost exclusively by effi- 
cient charged-current reactions in the accretion region perturbed 
by non-radial fluid flows. The spectrogram of the two phases 
is characterized by initially very small-scale angular variations 
with Z <: 12, which are associated with the onset of the Rayleigh- 
Taylor overturn activity, and which merge to continuously larger 
angular structures that correspond to I ^ 1 ... 4 modes toward 
the end of the accretion period at 0.4 - 0.6 s (depending on the 
model). This evolution is accompanied by a steady decrease of 
C™"" to a level of ~ 10"^ and a reduction of the electron neutrino 
flux density variations from values well beyond 100% to a level 
of ~ 10%, only (see Fig.|6] left panels). 



When neutrino-energy deposition in the post-shock lay- 
ers becomes sufficiently strong and the explosion is eventually 
launched at about 250 to 500 msec (depending on the model; 
Table 1), subsequent radial shock expansion rapidly diminishes 
the activity of the SASl and freezes post-shock convection. 
Single, longer lasting downdrafts of accretion flows are asso- 
ciated with isolated hot spots, where the variations of the total 
flux density can reach peak amplitudes of up to ~ 70 % (Fig.|5] 
panel 3). When accretion has ended, the amplitude of the angu- 
lar variations of the total neutrino energy flux reduces to a level 
of a few percent (Figs.[3]|4]l, and the angular pattern of the emis- 
sion becomes more uniform over the sphere, consisting of many 
spots with an angular size of ~ 30° (Fig.|5] panel 4). 
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Fig. 7. Pseudo-power coefficients C"^" of the electron-neutrino 
flux density as functions of angular mode number / at 200 ms 
(blue), 400 ms (red), and 1000 ms (black) for models W15-4 
(top), L15-3 (middle), and N20-3 (bottom), respectively. 

In the early post-accretion phase of model W15-4, 0.6 s 
t ^ 0.8 s, the spectrogram indicates the presence of low- 
amphtude (C"™ k lO"'*), small-angular size (Z ^ 10) pertur- 
bations in the electron neutrino energy flux caused by some low- 
amplitude turbulent flow in and below the neutrinospheric re- 
gion. When strong convection inside the PNS is encountered for 
t ^ 0.8 s the spectrogram drastically changes, being dominated 
by angular modes with I - 4, but stiU with C™'' ~ lO"'*. The 
electron neutrino flux density variations rise somewhat to a level 
of 10% to 20%, and become manifest in the total energy loss 
rate, too (Figs.[3]|4]l. 

Model N20-2 exhibits quite a similar behavior as model 
W15-4 except for the appearance of even larger (Z ~ 3) angular 
structures clearly recognizable in the pseudo-power spectrogram 
between 1.0 s and 1.2 s (Figs.|6ll2|. This differs from the behav- 
ior of model L15-3, where the amplitudes and angular size of the 
energy flux density variations remain small and even decrease in 
the post-explosion phase (Figs.|4l|6]l. 

The reason for the fluctuation behavior of the neutrino emis- 
sion during the vigorous pre-explosion and post-explosion ac- 




time [s] 

Fig. 8. Evolution of the non-radial specific kinetic energy (Vg + 
v^)/2 volume averaged over the computational domain inside 
the neutrinosphere for models W15-4 (solid), L15-3 (dashed), 
and N20-2 (dashed-dotted), respectively. 

cretion phases has been discussed, but what causes the spa- 
tial and temporal variations during the post-accretion phase? 
Because the explosion is well on its way at this time, neither 
post-shock convection nor the SASl nor accretion can be re- 
sponsible. Hence, there only remains non-radial gas flow in the 
outer layers of the proto-neutron star. Ledoux convection in the 
proto-neutron star thus may become visible eventually, i.e., its 
presence in the inner parts of the computational domain may 
become dominant in observable signals. This happens in mod- 
els W15-4 and N20-2, where the level of the non-radial specific 
kinetic energy (Vg + v^)/2, volume-averaged over the computa- 
tional domain, inside the neutrinosphere shows a steep rise at 
~ 0.8 s and ~ 0.9 s, respectively (Fig.[8]l. These non-radial flows 
that develop in models W15-4 and N20-2 at late times also be- 
come manifest in all discussed quantities: A^^^iQ., t), Amin(i^, t), 
C™", the dominant low /-modes {2^1% 4), and relative angu- 
lar flux-density variations. In contrast, no such eff'ect is present 
in model L15-3 (see Fig.|8]l, where we find a steady decrease 
of Amax(^^, - A„iin(il, 0. higher Z-modes {I ^ 10), smaller 
C™^", and lower flux-density variation amplitudes than in models 
W15-4 and N20-2 (see Figs.SS and©. 

Simulations with fully self-consistent treatment of the PNS 
interior show the presence of c onvection inside the PNS, 
i.e., b e low the neutrinosph ere (see iKeil et all 119961 : iBuras et al] 
l2006t iDessart et all l2006h more or less from the early post- 
bounce phase on. With the use of our inner radial grid boundary 
excising the inner parts of the PNS, and imposing neutrino lu- 
minosities at this boundary, convective activity is triggered only 
when the neutrino energy (or lepton number) inflow into the lay- 
ers close to the grid boundary is faster than neutrino transport 
can carry away this energy (or lepton number). Then convec- 
tively unstable gradients develop and convective flows begin to 
carry the energy and lepton-number outward. Whether this hap- 
pens or not depends on the boundary luminosities as well as on 
the location of the grid boundary within the density and tem- 
perature profiles of the PNS layers below the neutrinosphere. 
That location determines the efficiency of the neutrino transport 
and varies with the stellar progenitor, whose mass-infall rate de- 
cides how much mass accumulates in the near-surface layers of 
the PNS outside the inner grid boundary. The relative strength 
of the artificially imposed inflow of neutrino energy and lepton 
number compared to the efficiency of the neutrino transport on 
the grid, both sensitive to the location and contraction of the grid 
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electron neutrinos electron anti-neutrinos 




t [sec] t [sec] 

Fig. 9. Observable luminosity L,, (top row), observable mean en- 
ergy {E)o (middle row), and normalized quantity Lo{E}l (bot- 
tom row) of electron neutrinos (left column) and electron anti- 
neutrinos (right column) as a function of time for three of our 
models. Although we only present the results for one particular 
observer direction here, the global behavior and characteristics 
are very similar for all viewing directions. 



visible to an observer (Eq.fTOt. In Fig.|9]we show the observ- 
able electron neutrino and electron anti-neutrino luminosities for 
one chosen viewing direction for the three models W15-4, L15- 
3, and N20-2, respectively. The results for other directions look 
very similar with all characteristic features being independent of 
the observer position. We provide these quantities in addition to 
the total neutrino energy loss rate (Eqs.[T3] andfTSl Figs. [3] and 
lU, because their temporal evolutions are the ones expected to 
be measurable in the IceCube and Super-Kamiokande detectors. 
These detectors (mainly for Ve) will be sensitive to a combina- 
tion of the observable neutrino luminosity Lg and the observable 
mean neutrino energy {E}o- Thus, we also provide in Fig.|9]the 
time evolution of the observable mean neutrino energy and of 
the combination Lo{E}l, which (roughly) enters the IceCube de- 
tection rate of Cherenkov photons originating from the dominan t 
inverse beta decay reaction Ve + p ^ n + (iLund et alJIZOToh . 
Q Again one can recognize the different evolution stages, and 
in particular the post-shock convection and SASI phase, during 
which the quantity Lo{E}l exhibits rapid low-amplitude varia- 
tions for all three models. The level of the variations is a few per- 
cent (Fig.|9]l, which is considerably lower than that of the angular 
fluctuation amplitudes of the flux density, which reaches almost 
100% for the total neutrino flux density (Fig.|5]) and almost 200% 
for the electron neutrino and electron antineutrino flux densities 
(Fig-ISl lower panels). However, because the flux density vari- 
ations are caused by a few individual hot spots covering only 
angular areas of size ~ {nl9)^, the observable fluctuations (of Lq 
and {E)q) are lower by a factor of roughly (;r/9)^/(27r) ~ 1 /50. 
Some of this activity is also present at late times in the two mod- 
els W15-4 and N20-2, where Ledoux convection develops in the 
simulated outer parts of the proto-neutron star (see discussion 
above). 

From the results presented above we conclude that the sig- 
nals carry clear information about the postshock hydrodynamic 
activity, and about the duration and decay of the accretion pe- 
riod. Composition-shell interfaces present in the progenitor star 
can also leave an imprint. In model W15-4 the transition from 
the Fe-core to the Si-shell manifests itself in fast drops of the lu- 
minosities of Ve and at ~ 150 msec, when the mass accretion 
rate decreases steeply at the time the interface between the Fe- 
core and the Si-shell of the ISM© progenitor falls through the 
shock. 



boundary on the one hand and the chosen values of the boundary 
luminosities on the other, therefore decides about when, where, 
and how strongly convective activity develops below the neutri- 
nosphere. 

Because the position of and the conditions imposed at the 
inner boundary can thus influence the neutrino emission prop- 
erties, in particular during the post-accretion phase, our respec- 
tive model predictions must be considered with care. While they 
do not allow us to make any definite statements concerning the 
neutrino signal of a particular progenitor model because of the 
neglected treatment of the inner parts of the proto-neutron star, 
the models nevertheless show that convective flows below the 
neutrinosphere are likely to imprint themselves on the neutrino 
emission, and hence also on the GW signal of core-collapse su- 
pernovae. A measurement of these signals may actually provide 
some insight into the conditions inside proto-neutron stars. 

Because the neutrino energy flux density varies in our mod- 
els both with latitude and longitude, the observable neutrino lu- 
minosity Lo{t) is obtained by an integral over the hemisphere 



4. Gravitational wave signature 

Non-radial mass motions caused by gravity waves in the near- 
surface layers of the PNS, which are caused by the SASI 
and convection in the post-shock region as well as by convec- 
tive activity insid e the proto-neutron star (Murphv et al. 200^ 
iMarek et alj |2009) (driven by Ledoux unstable lepton or entropy 
gradients) result in a time-dependent, aspherical density strat- 
ification that produces gravitational radiation. The anisotropic 
emission of neutrinos associated with the non-radial mass flow 
(see Sect.O contributes to the gravitational wave signal, too. We 
computed and analyzed the signature of this gravitational radia- 
tion for the 3D models discussed in Sect. 12. 31 



"* Note that our transport approximation only provides luminosities 
and mean energies, but not the higher moments of the energy spectrum 
(see Sect.[3nj- 
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4.1. Formalism 

4.1.1. Non-radial mass flow 

If a source is of genuine three-dimensional nature, as it is the 
case for our models, it is common to express the gravitational 
quadrupole radiation tensor, h^, in the transverse traceless 
gauge in the following tensorial form 



1 



h''(X,t)^-(A^e^+A^e^) 

K 



(19) 



(see, e.g., iMisner et al. I [1991 . R denotes the distance between 
the observer and the source, and the unit linear-polarization ten- 
sors are given by 



(20) 
(21) 



with ee and being the unit polarization vectors in and 0- 
direction of a spherical coordinate system, and ® denoting the 
tensor product. 

The wave amplitudes A+ and Ax represent the only two in- 
depen dent modes of polarization in the TT gauge ( Misner et al.l 
Il993 h. In the slow-motion limit, they are obtained from lin- 
ear combinations of the second time derivatives (evaluated at 
retarded time, and denoted by a double dot accent) of the 
components of the transverse traceless mass quadrupole tensor 
(iMisner et al ] ll993h 



A+ = Qee-i 



(22) 
(23) 



We computed the latter using a post-Newtonian approach 
whereby the numerically troublesome second-order time deriva- 
tives of the mass quadrupole tensor components are trans- 
formed into much be tt er trac table spatial derivatives. Following 
iNakam ura & Ooharal (Il989h and iBlanchet. Damour & Schaferl 
(1990), the second-order time derivatives read in a Cartesian or- 
thonormal basis (the spatial indices / and j run from 1 to 3) 



7 = ^ J A\p [2viVj - Xi dj^^ff - xj 5,Oeif) , (24) 



where G is Newton's gravitational constant, c the speed of light 
in vacuum, the effective Newtonian gravitational poten- 
tial including the general relativistic "case A" correction of the 
monopole term according to Mareketal. (2006), p the mass- 
density, V, the Cartesian velocity components, and 5, the partial 
derivative with respect to the coordinate x' of a Cartesian basis. 

We note that the integrand in Eq. ( l24l ) has compact support 
and is known to the (2nd order) accuracy level of the numerical 
scheme employed in the hydrodynamics code. It can easily be 
shown that evaluating the integral of Eq. (|24|) by an integration 
scheme (of at least 2nd order) is by one order of accuracy supe- 
rior to twice applying numerical time-differentiati on methods to 
quadrupole data given at discre te points of time dPinn & Evan si 
ll990HMonchmeveret alJll99lh . 

Exploiting the coordinate transformation between the or- 
thonormal Cartesian basis x' and the orthonormal basis in spher- 
ical coordinates x' (with x' e [r, 0, 0]), the wave amplitudes A+ 
and Ax (Eqs. ( l22b & ( l23l l) are obtained from the following sec- 
ond time derivatives of the spherical components of the mass 



quadrupole tensor (lOohara et al.lll997tlScheidegger et al.ll2008l) 



^^COs2, 



• + ill sin^ ( 



+ 2 1^^ sin d> cos ( 



+l" sin^ 0-2 {llj cos (j) + I'^J sin (j)) sin 6* cos 6* , (25) 
ill ^ /I/sin^0 + /,^^cos^0-2/J,^sin0cos0, (26) 
llj = [l^J - /JJ) cos 6 sin (f>cos(f> + I^J cos 9 (cos^ 

- sin^ 0) + l].I sin sin - I^]^ sin 6 cos (p , 
where we used the abbreviation 



(27) 



TT 



(28) 



Cho osing (/> = one obtains the polarization modes (see, 
e.g.. lMisneretal.l (Il993h ) 



2^^ 



Ax 

for 9 -Q, and 



A, = lU-I^J 



(29) 
(30) 

(31) 
(32) 



for 9 - 7r/2, respectively. These expressions were already 
discussed in earlier investigations concerned with the evalua- 
tion of the gravitati onal wave signature o f 3D core-collapse 
supernova models (Muller & Jankal Il997t iFrver et all l2004t 
[Scheidegger et al. 2008, 201^! 

The total energy radiated in the form of gravitational waves 
due to nonspherical mass flow is given in the quadrupole approx- 
imation by (see, e.g., IMisner et al.ll 19931) 



C I ^ fTT ^xX^tTT 

„ "'[<''' 

./tTT i TT^2 . h i TT^2 . ,j TT^2 . 
+ (At ~ hz ) + " [y^xy ) + (4; ) + 



15G , 



jTT^^(jrT_jTT^ 



(33) 

(AD^)] , 



rr _ Q r TT 

ij ~ ij 

density is given by (where v denotes the frequency) 



with I.j = dl/j /dt, and the corresponding GW spectral energy 



dEu 2c^ 2 
^ ^ 15G^"""^ 



ITTT _ Trrp iTTT _ TTT]^ 



(34) 



+ \TV - Trt + 6 ( Wjf + \TI^f + \Ti 



where 



00 



^^(f)e"^'"'"df 



(35) 



is the Fourier transform of lf,^(t). 



4.1.2. Anisotropic neutrino emission 

To determine the gravitational wave signal ass ociated with the 
anisotropic emission of n eutrinos, we follow iMiiller & Jankal 
(Il997h and use Eq. (16) of lEpsteinI (Il978h in the hmit of a dis- 
tant source, /? — > 00, together with the approximation that the 
gravitational wave signal measured by an observer at time t is 
caused only by radiation emitted at time f' - t - R/c. Hence, we 
take t - f - const = Rjc, i.e., we assume that only the neutrino 
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Fig. 10. Relation between the source coordinate system 
(x' ,y' ,z') and the observer coordinate system {x,y,z). Changing 
from the observer system to the source system involves a rota- 
tion by an angle a about the z'-axis to an intermediate coordinate 
system (x*, y*, z*), followed by a rotation by an angle /? about the 
y*-axis (which thus is also the y-axis). 



pulse itself causes a gravitational wave signal, whereas memory 
effects, which prevail after the pulse has passed the observer, are 
disregarded. 

With these simplifications, the dimensionless gravitational 
wave amplitudes of the two polarization modes are given in the 
transverse-traceless gauge for an observer located at a distance R 
along the z-axis of the observer frame bV iMiiller & Janka (.19971) 

h+{t)^— \ dt' I dQ'(l -I- cos 61) cos 20 (fi',?') (36) 

C^R Jo J At, dii' 



and 



= ^ df' dQ'(l+cos0)sin2^*— (n',f'), (37) 

C^R Jo J4n 



respectively. Here dA(Q', f')/dn' is given by Eg. (fT2l l and de- 
notes the total neutrino energy radiated at time t' per unit of time 
into a solid angle dO' in direction (ff, (p'). Except for position- 
dependent factors the gravitational wave amplitudes are simply 
a function of this quantity provided by the ray-by-ray transport 
approximation (note that in Miiller & Janka (1997) we used the 
symbol Ly instead of A). 

The angular integration, dO.' = -d{cosO')d<p', in Eqs. (|36]| 
and ( |37] ) extends over all angles ff and <p' in the coordinate 
frame of the source {x',y',z') that we identify with the (arbitrar- 
ily chosen) spherical polar coordinate frame to which the hydro- 
dynamic results were mapped from the Yin- Yang grid employed 
in the simulations. For the evaluation of the polarisation modes 
we used the (asymptotic) values of dA(Q, f)/dQ extracted at a 
radius of 500km from our 3D models. 

The angles and in Eqs. ( |36] | and ( |37] ) are measured in the 
observer frame (x, y, z), while the neutrino luminosity is mea- 
sured in the source frame (x',y',z')- To allow for an arbitrary 
orientation of the observer relative to the source, we introduce 
two viewing angles a € [-7t,+7t] and /3 e [0,7r] (see Fig.fTOb. 
The coordinates measured in the observer frame are then related 
to the coordinates in the source frame by the following coordi- 
nate transformations 



X = X cos a + y sm a , 
y* — —x' sin a +y' cos a , 
z = z , 



(38) 
(39) 
(40) 



and 

X - X* cosyS - z* sinyS, (41) 

y = y* , (42) 

z = X* sinyS -I- z* cosy8. (43) 

With these coordinate transformations and the relations 

X - rsinflcos^, (44) 

y - rsin0sin0, (45) 

z = rcos0 (46) 

between Cartesian coordinates (x, y, z) and spherical polar coor- 
dinates (r, Q, 0), we obtain 

sin Q cos <p - (cos 0' cos a + sin <p' sin a) sin cos ji - 

cos 0' sin/?, (47) 

sin ^? sin - (sin <p' cos a - cos <p' sin a) sin & , (48) 
cos Q - (cos 0' cos a + sin <p' sin a) sin 6*' sin -i- 

cose'coSjS. (49) 

These expressions relate the angular coordinates in the observer 
frame (Q, 0) to those in the source frame { , (p'). For the specia l 
case a = they were already presented bv lKotake et alJ (l2009ah . 
Using Eq. ( |49] l and the equalities 



sin 2(p 



cos 20 = 



2xy 

+ y"^ 

2 2 
X - y 

x} + y"^ 



(50) 



(51) 



derived from Eqs. (l44l) to ( l46b . the two polarization modes 
(Eas.[36land[37]) are given by 

h{t,a,li)^— df'A(f')as(f',«,j8), (52) 
c^R Jo 

where S e (+, x) and A(f) is the angular integral of the neutrino 
energy radiated at time t per unit of time given in Eqs. ( fT3l l and 

as(f,a,/?)= r dQ'Ws(Q',a,/?)^(n',f), (53) 

are anisotropy parameters, which provide a quantitative measure 
of the time-dependent anisotropy of the emission in both polar- 
ization modes. Note that the evaluatio n of the ani sotropy param- 
eter ait) defined in Eq. (29) of Mull er & Jankal ([T997,) . which 
should not be confused with the observer angle a introduced in 
Fig. [To] does neither involve a dependence on observer angles 
{a,IJ) nor on the polarization mode. 

The angular weight functions appearing in the above expres- 
sion for the anisotropy parameters are given by 

Ws(6 ,(p ,a,p) - 



N{e',(l>',a,/3) 
where 

D+ = [I + (cos <p' cos a + sin (p' sin a) sin 0' sin /3 + cos 9' 
coSjS] {[(cos <f)' cos a + sin cp' sin a) sin 6' cos/5- 
cos 6' sin j8]^ - sin^ 6' (sin (p' cos a - cos (p' sin c)^ j 

Dx = [I + (cos <p' cos a + sin <p' sin a) sin 6' sin /3 + cos 6' 
cos p] 2 [(cos (p' cos a + sin cp' sin a) sin 6' cos /3 - 
cos 6' sin 0] sin 6' (sin (p' cos a - cos (p' sin a) , 
= [(cos 0' cos Q- + sin 0' sin Q-) sin 0' cosyS - cos 0' 
sin 0]^ + sin^ 6' (sin (p' cos a - cos (p' sin of . 



(54) 



,(55) 

(56) 
(57) 
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Choosing a - and j3 - 7t/2 the observer is located in the 
equatorial plane of the source (i.e., perpendicular to the source's 
z'-axis) at the azimuthal position 0' = 0. In that case one ob- 
tains simpler expressions for the angular functions (see also 
iKotake et aljr2009a) 



9 , , , , 1 + sin 6*' cos (h' 



W+l, = (cos^ 0' - sin-" 0' sin^ 0) — 



VKv L = -2 cos sin 6' sin S' 



cos- 0' + sin^ 9' sin^ (j)' 
1 + sin 9' cos (p' 
cos^ 9' + sin^ 9' sin^ ch' 



(59) 



Note that for axisymmetric sources hy, - 0. 

In general, the total energy Ecwit) radiated to infinity by a 
source in form of gravitatio nal waves until time t is given by 
(see, e.g. jMisner et alj (11993 ): Greek indices run from to 3, 
and repeated indices are summed over) 



£Gw(f)= r df' r TovnV^dQ, 
Jo Jsl 



(60) 



where the angular integration is performed over a two-sphere at 
spatial infinity 5^, and n'' = (0, 1, 0, 0) is a unit spacelike vector 
in polar coordinates {ct, r, 9, (p) normal to 5'^. Denoting by (. . . ) 
an average over several wavelengths, the gravitational-wave 
energy-momentum tensor r^y is given in transverse-traceless 
gauge by 



~ 32nG 

Thus, Eq. ( l60b can be rewritten as 



(61) 



Ecw(t) = ^ j| dt' r^dn{(d,hY)(drhi^)) , (62) 

where we have used the facts that /zJJ = 0, hJJ = 0, and 

c dyhJi^ = ~^rhJi^ for radially outgoing gravitational radiation. 
Evaluating the double sum in Eq. (l62b and using the r elations 
h^ff = -h-lJ = h+ and hj^J = /iJJ = hy, (see, e.g., .Misner et alJ 
Il993h . we finally find 

Ecwit) = dt' r^dQ {(drh^f + (d,hyf) . (63) 

Inserting the expressions for h+ and /ix given in Eq. ( |52] ) into 
Eq. (1631 ). we obtain for the energy £'n(0 radiated in form of grav- 
itational waves until time t due to anisotropic neutrino emission 

Enit) = ^ j| dr' dQ„^ ^1(1', a,/3) + ll(t', a,/3)] (64) 

with dQafi = sinySdySdor and 

kit,a,/3) = Ait)as{t,a,/3). (65) 
The corresponding spectral energy density is given by 



d£N G ij. ,2 
where l(v) is the Fourier transform of 



(66) 



/(f) = 1^ dQ„;, [ll(t, a,P) + ll(t, ff,^)]| . (67) 

For completeness we also provide an expression for the to- 
tal energy radiated in form of gravitational waves until time f. 



i.e., due to anisotropic mass flow and neutrino emission. It is 
obtained by inserting the total GW amplitude, i.e., the sum of 
the amplitudes given by Eqs. ( fT9] l and ( l52b into Eq. (|63] ). which 
leads to 



16;rG 



Jo J4n 



dD. 



.4 'x 



(68) 



4.2. Results 



Although an observer can only measure the total gravitational 
wave amplitude, i.e., that due to the combined effect of non- 
radial flow and anisotropic neutrino emission, we will first dis- 
cuss the GW signal of non-radial mass flow only, because it re- 
flects the various phases of the post-bounce evolution already 
introduced in the discussion of the neutrino signal above. 

Until post-shock convection and the S ASI are eventually ma- 
ture at around 150 msec, the GW signal is very small (Fig.fTTTi. 
Note that our models are not able to follow the GW emission 
that is caused by prompt post-sh ock convection bec ause of the 
excised inner region of the PNS dMarek et al.ll2009h . Later on, 
sizable g-mode activity is instigated in the outer layers of the 
proto-neutron star by convective overturn and the SASI dur- 
ing the hydrodynamically vigorous pre-explosion phase, and 
by the impact of anisotropic accretio n flows during the subse- 
quent post-explosion accretion phase (iMar ek et al.ll 2009 |) . This 
g-mode activity is the cause of GW sign als ( Marek etal]|2009l : 
iMurphv et all2009l: lYakunin et al.ll20i a). whose maximum am- 
plitudes are on the order of a few centimeters centered around 
zero. 

The GW frequency distribution possesses a very broad max- 
imum in the range of 100 Hz to 500 Hz, and the frequency corre- 
sponding to this maximum slowly increases with time (Fig.fTTTi. 
Partially already during the post-explosion accretion phase, but 
at latest when the shock wave starts to rapidly propagate to large 
radii between ~ 0.4 sec and ~ 0.7 sec (see Figs.|3]and|4]i, the GW 
amplitudes start to grow by about a factor of ten until approxi- 
mately asymptoting at ~ 0.9 sec in the case of the models based 
on the progenitor W 15, and at ~ 1 sec in the case of models based 
on the progenitor LI 5, respectively (Fig.fTTTi. This growth of the 
amplitude is associated with the anisotropic expansion of the 
shock wave, and a positive/negative w ave amplitude indic ates 
a prolate/oblate explosion, respectively (iMurphv et alJf2009t) . 

While the GW amplitudes grow, the GW energy distribu- 
tion dZiM/dv becomes narrower and dimmer, and the frequency 
at maximum power continues to increase. The latter effect was 
also observed in the 2D models of Murphy et al., (2009). At late 
times, the GW signal of the W15 models clearly signifies the 
convective activity inside the proto-neutron star through low- 
amplitude, high-frequency fluctuations around the asymptoti- 
cally roughly constant mean GW amplitudes, while no such fluc- 
tuations are present in the case of the L15 models (see discussion 
of the neutrino signal above). This model discrepancy is also 
evident from the energy spectrograms, which do exhibit a pro- 
nounced broad maximum (between ~ 350 Hz and ~ 550 Hz) 
at f > 0.8 sec in the case of the W15 models, but none for the 
L15 ones. Furthermore note that until the end of the simulations 
the frequency of the maximum of d/iM/dv has increased from 
around 100 Hz to almost 500 Hz for the former models (owing 
to the increasing speeds of mass motions in the postshock region 
at times $ 0.5 sec and because of the increased compactness of 
the proto-neutron star at times ^ 0.6 sec, respectively). 
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Fig. 11. The four panels show the gravitational wave amplitudes (top) and spectrograms of d^M/dv (bottom; normalized to the 
absolute maximum) arising from non-spherical mass flow of models W15-2 (top left), W15-4 (top right), L15-2 (bottom left), and 
L15-3 (bottom right), respectively. Blue curves give the amplitude A+ at the pole (solid) and the equator (dotted), while red curves 
show the other independent mode of polarization Ax from the same directions. 



The behavior of the total (matter plus neutrinos) GW am- 
plitudes is significantly different from that of the flow-only GW 
amplitudes for models that exhibit PNS convection below the 
neutrinosphere, i.e., for the models based on the progenitors 
W15 and N20. Particularly at late times, anisotropic neutrino 
emission causes a continuing growth of the GW amplitudes (in- 
stead of a saturation) in these models, while this is not the case 
for the L15 models (see Fig. [12] and compare with Fig.fTTTi. The 
latter behavior is also reflected in the time evolution of the asym- 
metry parameter as (Eq.l53Tl of the neutrino emission (Fig.[T3ll. 
The asymmetry parameter is practically zero in model L15-3 at 



late times, while it remains, after having temporarily grown to 
values beyond about 0.4 - 0.5%, at the level of ~ 0.3 % until the 
end of the simulations in models W15-4 and N20-2. 

The final GW amplitudes are up to a factor of two to three 
higher when taking the contribution of anisotropic neutrino 
emission into account. In contrast, the amount of energy radi- 
ated in the form of GW, which is proportional to the third time- 
derivative of the quadrupole moment and hence proportional to 
the time derivative of the GW amplitude is only insignificantly 
changed, and is practically constant for all simulated models af- 
ter the onset of the explosion (see Fig.[T4li. The integral value 
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LI 5-3 (flow + i^'s) 
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Fig. 12. Gravitational wave amplitudes Rh^^ (blue) and Rhl.^ (red) due to anisotropic mass flow and neutrino emission as a function 
of time for models W15-4 (left), L15-3 (middle), and N20-3 (right) , respectively. The solid curves show the amplitudes for an 
observer located above the north pole (a = /? = 0; see Fig.fTOli of the source, while the other curves give the amplitudes at the 
equator (a = 0, y6 = 7r/2). 
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Fig. 13. Asymmetry parameter as of the neutrino emission (Eq.l53il as a function of time for models W15-4 (left), L15-3 (middle), 
and N20-2 (right), respectively. The panels in the upper row show a+ (blue) and (red) for a particular observer direction, while 
the panels in the lower row give for both parameters the maximum and minimum values in all directions. 



of the GW energy radiated by neutrinos is low 1%) com- 
pared to that emitted by matter through the slow variation of 
the GW neutrino amplitude with time, i.e., its time derivative is 
much smaller than that of the GW matter amplitude. For this 
reason we also abstained from evaluating the total energy radi- 
ated in form of GW (Eq.l68Tl. It differs little from that caused by 
anisotropic matter flow alone (Eq.[33Tl. because the mixed term 
in Eq. ( |68] |. resulting from the square of the sum of the matter 
and neutrino parts, contributes ^ 10% to the total radiated GW 
energy, and the pure neutrino term ^ 1%. Figure [14] also shows 
that the (small) contribution of anisotropic neutrino emission to 



the radiated GW energy is enhanced at late times when proto- 
neutron star convection occurs below the neutrinosphere, as it is 
the case for models W15-4 and in particular N20-2. 

The variation of the total GW amplitudes with observer an- 
gle is illustrated in Fig.[T5]for model W15-4 at 1.3 sec (when the 
simulation was stopped). Both the amplitude variations and the 
typical angular size of the speckled GW emission are similar for 
all other simulated models. The model-independent level of the 
amplitude variations is also supported by Fig. [12] when compar- 
ing various amplitudes at any given (late) time. 
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Fig. 15. Gravitational wave amplitudes due to anisotropic mass 
flow and neutrino emission, Rh^^ (top) and Rh^^ (bottom), as 
functions of the observer angles (see Fig.fTOll for model W15- 
4 at 1.3 sec past bounce. The white contours give the locations, 
where the amplitudes are zero. Yellow and red areas indicate 
positive amplitudes, green and blue negative ones. 



Fig. 14. Energy emitted in form of gravitational waves due to 
anisotropic mass flow (top panel) and due to anisotropic neutrino 
emission (bottom panel) as a function of time for models W15-4 
(solid), L15-3 (dashed), and N20-2 (dash-dotted), respectively. 



The (normalized) amplitude spectrograms of the total grav- 
itational wave amplitudes (d(A+.x + Rh+,x)/dv; Figs.fTSlandfTTTi 
illustrate two model-independent findings. Firstly, during the hy- 
drodynamically vigorous pre-explosion and post-explosion ac- 
cretion phases (0.2 ^ t $ 0.5 - 0.7 sec) the spectra of all models 
are characterized by some power at low frequencies (^ 100 Hz) 
and a broad power maximum at frequencies ~ 200 Hz and an- 
other weak one at ~ 800 Hz. The latter broad maximum at 
high frequency is more pronounced in the models based on the 
W15 and N20 progenitors and in the cross polarization GW 
mode. Secondly, during the post-accretion phase {t ^ 0.7 sec) 
the spectra of all models are dominated by a low-frequency 
(^ 40 Hz) contribution peaked toward the lower end of the spec- 
trogram. In the models where PNS convection occurs below the 
neutrinosphere (models W15 and N20) we also find a double- 
peaked high-frequency contribution decreasing/increasing from 
~ 700 Hz (400 Hz) at f ~ 0.8 sec, and eventually merging into a 
single power maximum at ~ 500 Hz at f ~ 1.2 sec. Again, this 
contribution is more pronounced for the cross polarization GW 
mode. 

The spectra of the total GW amplitudes are dominated by 
the contribution from non-isotropic neutrino emission at low 
(< 100 Hz) frequencies (Figs.fTSl and [TTli. At higher frequen- 
cies (> 100 Hz) the spectra of model W15-4 show two pro- 
nounced maxima (at 100 - 200 Hz and 600 - 800 Hz, respec- 
tively) at all times. These maxima are also present in model 



L15-3 at times < 0.7 s, the high-frequency one being, however, 
much less pronounced. The lower maximum (at 100 - 200 Hz) 
results from g-mode activity in the PNS surface instigated by 
non-radial flow (SASI, accretion) in the post-shock region until 
~ 0.5 - 0.7 sec. At later times PNS convection is responsible for 
the peak between 300 and 500 Hz. We have proposed this ex- 
planation already for the corresponding maxima present in the 
GW energy spectrograms arising from non-spherical mass flow 
(Fig.fTTTi. and discussed why the frequencies of these maxima 
increase with time. The source of the high-frequency maximum 
(600 - 800 Hz) is unclear, but a further detailed analysis shows 
that (i) the maximum is solely caused by non-radial gas flow, 
i.e., it is not connected to neutrinos, (ii) it does not result from 
stellar layers below the neutrino sphere but from those close to 
or slightly above it, and (iii) does not depend on the position of 
the observer 

Note that the high-frequency maximum present in the am- 
plitude spectrograms is strongly suppressed in the correspond- 
ing energy spectrograms (Fig.fTTTi. because the latter involve the 
squared time derivatives of the amplitudes. Thus, the already 
high ratio of the low- and high-frequency maxima in the am- 
plitude spectrograms (about two orders of magnitude) translates 
into an even higher ratio for the energy spectrogram maxima, 
rendering the high-frequency maximum practically invisible. 

5. Discussion and conclusions 

Based on a set of three-dimensional (3D) parametrized neutrino- 
driven supernova explosion models of non-rotating 15 and 
20 Mq stars, employing a neutrino transport description with 
a gray spectral treatment and a ray-by-ray approximation of 
multi-dimensional effects (the scheme is applicable in the regime 
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Fig. 16. Normalized (to the absolute maximum) amplitude spectrograms of the total gravitational wave amplitudes A+ + Rh+ (left 
panels) and Ax + Rh^ (right panels) at the pole for model W15-4. The lower panels show the spectrograms in the frequency range 
5 Hz to 100 Hz, and the upper ones in the frequency range 100 Hz to 1 kHz. 
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Fig. 17. Same as Fig.[T6l but for model LI 5-3. 



outside the dense neutron star core, i.e., around and outside 
the neutrinosphere), we evaluated both the time-dependent and 
direction-dependent neutrino and gravitational-wave emission of 
these models. To this end we presented the formalism necessary 
to compute both the observable neutrino and gravitational wave 
signals for a three-dimensional, spherical source. For the neu- 
trino signal we presented formulas that allow one to estimate the 
apparent luminosity when the local flux density F on a sphere is 
known and the ansatz of Eq. (4) about the angular distribution of 
the intensity can be made. While in general the location of the 
neutrino-decoupling surface has to be suitably defined {e.g., by 
a criterion based on the optical depth), our ray-by-ray transport 
scheme implies that r^F =const in the free-streaming limit in 
every direction. Thus, it spares us such a definition of the neu- 
trinosphere for evaluating the angular integration for the observ- 
able luminosity. Concerning the gravitational-wave analysis, we 
extended and generalized previous studies, where the source was 
either assumed to be axisymmetric or where the formulas for the 



signals of a 3D source were only given for special observer di- 
rections. 

Our models followed the evolution from shortly after core 
bounce up to more than one second into the early cooling evo- 
lution of the PNS without imposing any symmetry restrictions 
and covering a full sphere. The extension over such a relatively 
long evolution time in 3D was possible through the usage of an 
axis-free overset grid (the Yin- Yang grid) in spherical polar co- 
ordinates, which considerably eases the CFL time-step restric- 
tion and avoids axis artifacts. A central region, the dense inner 
core of the proto-neutron star, was excised from the computa- 
tional domain and replaced by an inner, time-dependent radial 
boundary condition and a gravitating point mass at the coordi- 
nate origin. Explosions in the models were initiated by neutrino 
heating at a rate that depends on suitably chosen values of the 
neutrino luminosities imposed at the inner radial boundary. 

The post-bounce evolution of our models can be divided into 
four distinct phases (Fig. [3]). The first phase, the quasi- spherical 



18 



E. Miiller et al.: Parametrized 3D models of neutrino-driven supernova explosions 



shock-expansion phase, lasts from shock formation shortly after 
core bounce to 80 - 150 msec, when convection sets in. The sec- 
ond phase, the hydrodynamically vigorous pre-explosion phase, 
comprises the growth of post-shock convection and of the stand- 
ing accretion shock instability (SASI). The post-explosion ac- 
cretion phase begins when energy deposition by v-heating in 
the post-shock layers becomes sufficiently strong so that the to- 
tal energy in the post-shock region ultimately becomes positive. 
During this phase the shock accelerates outward while gas is still 
accreted onto the PNS. This process is commonly called "shock 
revival". The duration of the latter two phases depends on the 
progenitor. During the post-accretion phase, the fourth and final 
phase characterizing the evolution of our models, accretion ends 
and the proto-neutron star develops a nearly spherical neutrino- 
driven wind. 

The neutrino emission properties (fluxes and effective spec- 
tral temperatures) of our 3D models exhibit the generic time- 
dependent features already known from 2D (ax isymmetric) 
models (e.g.. .Buras et al. 2006: Scheck et ani2006t [Marek et aL. 
l2009tlBrandt etalJl201 lb, showing fluctuations over the neutron 
star surface on different spatial and temporal scales. We found 
that non-radial mass motions caused by the SASI and convection 
in the neutrino-heated hot-bubble region as well as by PNS con- 
vection below the neutrinosphere give rise to a time-dependent, 
anisotropic emission of neutrinos, particularly of electron neu- 
trinos and anti-neutrinos, and thus also to the emission of grav- 
itational waves. Because very prominent, quasi-periodic slosh- 
ing motions of the shock due to the standing accretion-shock 
instability as visible in 2D simulations are absent and the emis- 
sion from different surface areas facing an observer adds up in- 
coherently, the modulation amplitudes of the measurable neu- 
trino luminosities and mean energies are sign ificantly lower than 
predicted by 2D m odels (for 2D results see iMarek et alj|2009l; 
iBrandt et alJl20TTh . 

During the quasi-spherical shock expansion phase shortly af- 
ter bounce the level of temporal and angular fluctuations of the 
neutrino emission is low ($ 10"^). In contrast, the fluctuation 
amplitudes reach a level of several 10% of the average values 
during the hydrodynamically vigorous pre-explosion phase and 
the post-explosion accretion phase, where a few distinct, highly 
time-variable regions or even short-lived single spots with an an- 
gular size of 10° to 20° are responsible for the brightest emission 
maxima. As the outward shock expansion is well on its way in 
the post-explosion accretion phase, still existing accretion down- 
drafts can be responsible for similar fluctuations in the neutrino 
emission, though the number of corresponding hot spots de- 
creases with diminishing accretion. When accretion has ended 
and the post-accretion phase has started, directional variations 
can be caused by the occurrence of Ledoux convection in the 
outer layers of the proto-neutron star, which we indeed observe 
in models based on two of our three progenitors (see also the 
discussion of the influence of the inner radial boundary condi- 
tion below). The temporal and angular variations of the emission 
in different directions are even more pronounced when consid- 
ering the energy flux of the electron neutrinos or electron anti- 
neutrinos alone (instead of the emission in all neutrino flavors). 
In that case the angular variations of local flux densities can 
exceed 100% in all models during the pre-explosion and post- 
explosion accretion phases, and the peak values can be close to 
200% during short episodes. The total energy loss rates in neutri- 
nos and the observable luminosities as surface-integrated quan- 
tities, however, are much smoother in time during all phases, 
showing fluctuation amplitudes of at most several percent. 



The gravitational wave emission also exhibits the generic 
time-dependent features already known from 2D (axisymmet- 
ric) models, but the 3D wave amplitudes are considerably 
low er (by a factor of 2 -3) than those p r edicted by 2D mod- 
els (Mu Uer et all 120041: iMarek et all 120091: iMurphv et aT]|2009l: 
YakuninetalJ|2010!) owing to less coherent mass motions and 
neutrino emission. Note in this respect that the GW quadrupole 
amplitudes, which are usually quoted for 2D models (Ajq ), have 
to be multiplied by a geometric factor sin^0 V15/7r/8 (which 
is equal to s; 0.27 for 9 - 90°)- Violent, non-radial hydrody- 
namic mass motions in the accretion layer and their interaction 
with the outer layers of the proto-neutron star give rise to a GW 
signal with an amplitude of ~ 2-4 cm in the frequency range 
of ~ 100 Hz to ~ 400 Hz, while anisotropic neutrino emission is 
responsible for a superimposed low-frequency evolution of the 
wave amplitude, which thus can grow to maximum values of 
10 - 20 cm. Variations of the mass-quadrupole moment caused 
by convective activity inside the nascent neutron star contribute a 
high-frequency component (300 - 600 Hz) to the GW signal dur- 
ing the post-accretion phase. The GW signals exhibit strong vari- 
ability between the two polarizations, different explosion simu- 
lations and different observer directions, and besides common 
basic features do not possess any template character. 

Finally we would like to reflect on some of the deficiencies 
of the presented 3D models. Because of the ray-by-ray treat- 
ment of the v-transport, the directional variations of the neutrino 
emission in respo nse to local inhomogeneities in th e star may 
be overestimated ( Ott et al.ll2008t iBrandt et al.ll201 ih . However, 
when we evaluate observable signals, these artificial effects are 
mostly compensated for by integrations of the neutrino flux den- 
sities over the surface areas visible to observers from different 
viewing directions (see Eqs.lSllTl fTOlllll i. or by the integration 
of the neutrino energy loss in all directions (Eqs. [T3ll53T l. 

Another deficiency concerns the usage of the inner radial 
grid boundary, because of which our simulations do not fully 
include (either in space or time) the convective flow occur ring in 
the PNS interior after co re bounce (iKeil et al 1996; Bura s et al] 
2006; Dessart et al.ll200 6). Moreover, convective activity in the 
simulated outer layers of the PNS occurs for special conditions: 
It is triggered only when the artificially imposed inflow of neu- 
trino energy and lepton number through the inner radial bound- 
ary into the adjacent layers is faster than the neutrinos can carry 
away this energy or lepton number. Whether this is the case 
sensitively depends on the employed neutrino-transport approx- 
imation, but also on the location and contraction of the grid 
boundary, the chosen values of the boundary luminosities, and 
on the stellar progenitor. Its mass-infall rate decides how much 
mass accumulates in the near-surface layers of the PNS outside 
the inner grid boundary. Because the position of and the condi- 
tions imposed at the inner boundary can thus influence the neu- 
trino emission properties, in particular during the post-accretion 
phase, our respective model predictions must be considered with 
care. While they do not allow us to make any definite statements 
concerning the detailed neutrino signal of a particular progeni- 
tor model due to the neglected treatment of the inner parts of the 
proto-neutron star, the models nevertheless show that convective 
flows below the neutrinosphere are likely to imprint themselves 
on the neutrino emission, and hence also on the GW signal of 
core-collapse supernovae. A measurement of these signals may 
actually provide some insight into the conditions inside proto- 
neutron stars. 
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